Astronomy & Astrophysics manuscript no. vojtech'simon-GRB090726-revl 


© ESQ 2009 


November 9, 2009 





Complicated variations of early optical afterglow of GRB 090726 

V. Simon 1 , C. Polasek 1 , M. Jelmek 2 , R. Hudec 1 , J. Strobl 1 

1 Astronomical Institute, Academy of Sciences of the Czech Republic, 251 65 Ondfejov, Czech Republic 

2 Instituto de Astrofisica de Andalucfa CSIC, Apartado de Correos, 3.004, E- 18.080 Granada, Spain 

Received date; accepted date 

ABSTRACT 

Aims. We report on a detection of an early rising phase of optical afterglow (OA) of a long GRB 090726. We resolve a complicated 
profile of the optical light curve. We also investigate the relation of the optical and X-ray emission of this event. 
Methods. We make use of the optical photometry of this OA obtained by the 0.5 m telescope of AI AS CR, supplemented by the data 
obtained by other observers, and the X-ray Swift/XKI data. 

Results. The optical emission peaked at ~ 17.5 mag(i?) at t - T » 500 s. We find a complex profile of the light curve during the 
early phase of this OA: an approximately power-law rise, a rapid transition to a plateau, a weak flare superimposed on the center 
of this plateau, and a slowly steepening early decline followed by a power-law decay. We discuss several possibilities to explain the 
short flare on the flat top of the optical light curve at t - To ~ 500 s; activity of the central engine is favored although reverse shock 
cannot be ruled out. We show that power-law outflow with o bs/®c > 2.5 is the best case for OA of GRB 090726. The initial Lorentz 
factor is To ~ 230 - 530 in case of propagation of the blast wave in a homogeneous medium, while propagation of this wave in a 
wind environment gives T ~ 80 - 300. The value of T in GRB 090726 thus falls into the lower half of the range observed in GRBs 
and it may even lie on the lower end. We also show that both the optical and X-ray emission decayed simultaneously and that the 
spectral profile from X-ray to the optical band did not vary. This is true for both the time period before and after the break in the X-ray 
light curve. This break can be regarded as achromatic. The available data show that neither the dust nor the gaseous component of the 
circumburst medium underwent any evolution during the decay of this OA, that is after t —T < 3000 s. We also show that this OA 
belongs to the least luminous ones in the phase of its power-law decay corresponding to that observed for the ensemble of OAs of 
long GRBs. 

Key words. Gamma rays: bursts - Radiation mechanisms: non-thermal - Plasmas - ISM: jets and outflows - Galaxies: ISM - 
Galaxies: starburst 



1. Introduction 

GRB 090726 was a long gamma-ray burst (GRB) localized by 
Swift/BAT on the 26 th July 2009 at 22:42:27.8 UT. It displayed a 
single peak profile, with Tgo (15-350 keV) = 67 ± 15 s. Fluence 
in the 15-150 keV band was (8.6+1.0)xl0~ 7 erg cirT 2 s _1 . The 
1-s peak flux was 0.7±0.2 photon crrT 2 s _1 (Page et al. 2009). 

The GRB coordinates were available via GCN within 65 s, 
although because of an Earth-limb constraint Swift could not 
follow-up until ~3.6ks after the trigger. Redshift of this GRB 
is z — 2.71, as determined from the lines in optical spectrum of 
the OA (Fatkhullin et al. 2009). 

In this letter, we report the early observations of optical 
afterlow (OA) of this GRB performed by the robotic 0.5 m 
telescope (D50) of the Astronomical Institute of AS CR in 
Ondfejov, Czech Republic. The first image of the D50 started at 
22:45:41 .0 UT, i.e. 194 s after the trigger. The telescope was tak- 
ing images until 23:35:35 UT (for 50 minutes), until a pointing 
limit was reached. On all of the 107 unfiltered, 20 s exposures 
the OA can be detected, although, to improve the signal to noise 
ratio, we co-added some of the images. 

2. Observations and data analysis 

D50 is a robotic, 500/1975 mm Newtonian telescope, located 
at AI AS CR in Ondfejov. It is equipped with a CCD cam- 
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Table 1. Calibration stars used, together with their USNO-A2.0 
derived magnitudes. Denoted error estimates (lcr) are for rela- 
tive brightnesses within this frame, external error - binding to 
USNO system (-0.04 mag) and error of the USNO calibration 
itself (~0.2 mag) - are much larger than this. 
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era FLI IMG 4710 (CCD chip E2V 47-10, mid-band coating, 
1024x1024 pixels), focuser FLI DF-2, and a field corrector 
TeleVue Paracorr PSB-1 1000. It is controlled by RTS2 software 
(Kubanek et al. 2006) and makes use of astrometric code by 
Lang et al. (2009). The field of view is 20 x 20 arcmin. In order 
to reach deeper magnitudes, unfiltered observations were carried 
out. Their peak sensitivity is close to the R filter. Exposure time 
of 20 s was used for each CCD image. 

Aperture photometry using phot routine within 
IRAF/DAOPhofl was used on the OA and 8 comparison 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
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Fig. 1. CCD image of OA of GRB 090726 secured with the 
D50 cm telescope. It was obtained close to the time of the peak 
brightness of this OA, at t - To » 670 s. Only a part of the orig- 
inal image is shown. OA and the comparison stars are marked. 
See Sect.|2]for details. 



stars. The reference grid of these 8 stars was then calibrated 
with brightnesses from USNO-A2.0 (Monet et al. 1998) using 
robust fitting algorithms to avoid possible influence of star 
variability. The values we obtained for brightnesses of the 
calibration stars are in Table [T] As we have no afterglow color 
information, we cannot transform the CR (clear, /^-calibrated) 
magnitudes into real R band. Because the airmass changed 
only between 1.180 and 1.124 during the observation, the 
atmospheric color effects should be negligible. CR magnitude of 
the OA differs from the R band only by the difference between 
the OA color and the mean color of the calibration stars. If 
we assume the (typically) constant color of the OA, also this 
difference stays constant. The data are given in Table [2] 

An example of our CCD image of the OA, obtained at t— Tq * 
670 s, that is near its peak magnitude, is displayed in Fig.[T] Only 
a part of the original 20x20 arcmin CCD image is shown. The 
OA is clearly defined as a new, variable source. 

The resulting light curve of the OA is displayed in Fig. [2] The 
R band observations of Moskvitin et al. (2009) and Volnova et al. 
(2009) were also included. Notice that the data by Moskvitin et 
al. (2009) fit our unfiltered observations nicely. The optical tran- 
sient was observed with a 60 cm telescope at Crni Vrh (Maticic 
& Skvarc 2009) simultaneously with our measurements. We 
found that their light curve shows a large scatter. The main con- 
tribution of their data is that their first point was obtained before 
the start of our observations (f - Tq = 154 s vs. t - To — 194 s). 
The brightness was fainter than or comparable to that of our first 
observation, which helps constrain any possible large-amplitude 
initial flash. 



for Research in Astronomy, Inc., under co-operative agreement with 
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Fig. 2. Light curve of OA of GRB 090726 in the red spectral 
region. The data coming from the individual observers are re- 
solved. HEC13 fit to Ondfejov data is marked by the smooth, 
solid line, (a) The very early phase plotted with a linear scale on 
the abscissa, (b) The early phase with a logarithmic scale of the 
abscissa, (c) Residuals of HEC13 fit. The open circles are high- 
lighting the flare, (d) The first day of OA. A vertical line denotes 
the start of SwiftPiKT observations. The error bars represent the 
uncertainties at the lcr level. See Sect.|2]for details. 
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Table 2. Observing log of Ondfejov data. Time elapsed from 
GRB trigger in seconds is marked as t - Tq while E refers 
to the exposure time in seconds (some images are co-added). 
Magnitude is denoted as mag, including its lcr error. 
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In order to resolve the profile of the early phase of OA in 
detail, both the panels with a linear and log-log scales of the 
axes are shown in Fig. [2] 

In order to show the profile of the early phase of OA more 
clearly and to lower the scatter of the observations, Ondfejov 
data were fitted by the code HEC 1 3 written by Harmanec (1992). 
The code is based on the method of Vondrak (1969, 1977), 
who improved the original method of Whittaker (Whittaker & 
Robinson 1946). Full description of the method can be found 
in Vondrak (1969). This method can fit a smooth curve to the 
data no matter what their profile is. These data can be distributed 
in time non-uniformly. The method is based on minimizing the 
value Q = F + A 2 S , where F — £ p(yi - y' { ) 2 denotes the de- 
gree of smoothing (y being the smoothed and y' the observed 
value of the variable), S = Yji^yd 2 is the measure of rough- 
ness of the curve, A 2 is a constant to be selected and defines how 
much the curve will be smoothed. HEC 13 makes use of two in- 
put parameters, e (in dimensionless units) and AT. The quantity 
e = 1 1 A 2 determines how "tight" the fit will be, that is if only 
the main profile or also the high-frequency variations are to be 
reproduced. The quantity AT is the interval over which the data 
are binned before smoothing. The resulting fit consists of the 
mean points, calculated to the individual observed points of the 
curve. A set of fits to the data with the different input parameters 
£ and AT was generated and submitted to inspection. Only the 
fit to the plateau and the decay is shown in Fig. [2] because the 
transition between the rising branch and the plateau was proved 
to be too sharp to give any reliable fit. Also the narrow flare at 
t - Tq as 500 s was omitted from the fitting. It was found that the 
fit with e=l, Ar = 400 s reproduces the main features of the 
light curve and significantly enhances the visibility of the main 
profile. The residuals of this fit are displayed in Fig.|2j; and show 
that four points, belonging to the flare and not included in the fit, 
clearly deviate from the profile of this curve. This fit also allows 
us to conclude that not other flare with the amplitude larger than 
the errors of the data occurred during Ondfejov observations. It 
is true that this method is somewhat subjective but it enables us 
to find a compromise between a curve running through all the 
observed values and an ideal smooth curve. We preferred to use 
this method because it does not make any assumptions about the 
profile of the fitted data. 
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Fig. 3. (a) Time evolution of the relation of the optical and X-ray 
flux (0.3-10 keV) in OA of GRB 090726, defined in the text, (b) 
Evolution of this relation with the X-ray flux, (c) Evolution of X- 
ray hardness H with the X-ray flux. The optical and X-ray flux 
of the OA were declining in all three diagrams. The error bars 
were calculated from the uncertainties given on both the optical 
and X-ray data at the 90 percent level. See Sect.[2]for details. 



The relation of the optical and X-ray flux, and the evolution 
of the X-ray hardness ratio were investigated using Swift/XRT 
(Evans et al. 2007, 2009) data of the afterglow of GRB 090726. 
Only the R band data of Moskvitin et al. (2009) were used be- 
cause Ondfejov observations must have already finished. In or- 
der to get an analogy to the color index, the quantity RX = 
log[F v (R)]-log[I x /v], where F V (R) is the flux in the R band, I x 
is the observed X-ray flux in the 0.3-10 keV band, and v is the 
frequency, was calculated. The conversion between the R band 
magnitude and flux was carried out according to the relation of 
Clark et al. (2000). The value of v = 4 x 10 17 Hz was used 
for Ix- Also hardness H of the X-ray spectrum was determined 
from the countrates cr in the hard and soft bands as follows: H = 
log(cr hard (1.5-10 keV))-log(cr soft (0.3-1.5 keV)). The results are 
displayed in Fig. [3] Notice that both RX and H remain constant 
within the observational errors during the decay of the afterglow. 

Spectral energy distribution (SED) of the afterglow of 
GRB 090726 is displayed in Fig. [4] The decay slope of Ondfejov 
data for t - Tq > 1730 s was extrapolated to t - Tq = 4000 s, 
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Fig. 4. SED of the afterglow of GRB 090726 at t - T = 4000 s. 
The X-ray data come from Swift/XRT, the optical measurement 
represents the extrapolation of Ondfejov data. See Sect.[2]for de- 
tails. 



that is to the time for which SwiftPiKT observations exist. The 
X-ray data were corrected for the absorption given by Page et 
al. (2009). The R band magnitude was corrected for the Galactic 
reddening by Ar — 0.1 1 mag according to the maps by Schlegel 
etal. (1998). 

3. Results 

We report on a detection of early, rising OA of GRB 090726. Its 
light curve displayed a complicated profile, which can be sum- 
marized as follows: 

1) Steep rise until t — Tq » 300 s which finished by a rapid 
transition to a plateau lasting until t — Tq » 800 s. 

2) A short flare on the flat top (plateau) of the optical light curve 
at t - T * 500 s, 

3) Only a slowly steepening decline in 800 s < t - Tq < 1400 s. 

4) The power-law decay steepened at t - Tq « 1400 s and lasted 
until at least till t - T » 7000 s. 

Another break in the decay occurred between t-To ~ 7000 s 
and t-To~ 0.86 d (see Volnova et al. 2009). 

The initial rise for t - To < 300 s (the first 4 points in 
Fig. [2} can be approximated by a power-law with the index 
a r - 0.94 + 0.05. This small a x speaks in favor of propagation 
of the blast wave in a wind of the progenitor rather than in the 
interstellar medium. Since the observed part of the rise of OA 
occurred well after the end of the gamma-ray emission, this is a 
thin shell expansion (e.g. Zhang & Kobayashi 2005). 

A short flare, superimposed on the flat top of the optical light 
curve, attained the peak magnitude at t - Tq w 500 s. A power- 
law fit yields a T p = 2.7±0.8 for the rising branch, while the decay 
can be fitted with = -3.6 ± 0.8. This decay is considerably 
steeper than the later power-law one. Prompt emission of GRB is 
ruled out because the flare occurred well after the end of the ob- 
served GRB emission. Also a density enhancement is less likely 
because the decline shallower than r^ 2+ ^' is expected (Fenimore 
et al. 1996). Assuming a common spectral index f3 = 1 (Nardini 
et al. 2006), adF should be < -3 in this case. If this flare is due 
to reverse shock, then it sorts OA of GRB 090726 to Type II, 
defined by Zhang et al. (2003). However, the observed decay ap- 
pears to be too steep because o^f ~ 2 is expected (Kobayashi 
2000), although a wind could help steepen this decline (Greiner 
et al. 2008). In this regard, we note that any additional flash, 
like the one in GRB 990123 (Akerlof et al. 1999) and caused by 
another reverse shock, could occur before the start of Ondfejov 
observations only if its peak magnitude were comparable with 
or fainter than that of the peak observed at t - To ~ 500 s, The 



absence of this flash is strengthened by the fact that also the first 
point of Maticic & Skvarc (2009) shows that the OA was not 
brighter than in the subsequent Ondfejov data. The steep decay 
of the observed flare speaks in favor of activity of the central 
engine (Fan et al. 2008). This flare occurred during the phase 
in which OAs observed in X-rays often display a plateau and/or 
flares (e.g. Burrows et al. 2005; Willingale et al. 2007). A flare 
can appear if the central engine continues to create shells of 
much smaller power at late times with smaller T even after GRB 
(Ghisellini et al. 2007). 

The profile of OA of GRB 090726 near the peak magnitude 
can be used to discriminate between various profiles of the jet 
and the viewing angles, using the recent models by Panaitescu & 
Vestrand (2008). The top of this OA is strikingly different from 
the sharp ones of GRB 060418 and GRB 060607 A (Molinari et 
al. 2007; Nysewander et al. 2009). Its flat top speaks against 
isotropic outflow or a jet seen face-on. Jet with a sharp angu- 
lar boundary with o bs/©jet > 2 can give rise to a short plateau, 
but the time of its peak luminosity f pea k occurs at t - To > 1000 s 
for z = 2.71, which is too late to be consistent with our data. 
Power-law outflow with o bs/©c > 2.5 (0 C being angular size of 
the core) is the best case for OA of GRB 090726. It yields a flat 
top at t - Tq of the order of several hundred seconds. The flat top 
thus can be caused by observing the power-law jet of the OA of 
GRB 090726 off-axis. 

Some OAs were observed to display a re-brightening or a 
plateau between rest frame (f - 7b) res t of 0.01 and 0.02 d (i.e. 
t - To between -3200 and -6400 s f or z = 2.71) (Klotz et al. 
2005). This time period corresponds to the phase of the power- 
law decay in OA of GRB 090726. Also the X-ray 0.3-10 keV 
data display only smooth, power-law decay without any plateau 
in this epoch. 

Determination of the initial Lorentz factor To of GRB re- 
quires a knowledge of its isotropic energy release E\ so . Since 
only a 15-150 keV spectrum, best fitted by a simple power- 
law model (Page et al. 2009), is available, we estimated E lso 
from Amati relation (Amati et al. 2002). The observed spec- 
trum suggests E[ so > 2 x 10 52 erg. The upper limit of GRBs 
is Eiso ~ 3 x 10 54 erg. Following the approach of Molinari et 
al. (2007), the initial Lorentz factor of GRB 090726 turns out to 
be To ~ 230 - 530 in case of propagation of the blast wave in a 
homogeneous medium, while propagation of this wave in a wind 
environment gives To ~ 80-300. The density n — 1 cirT 3 and the 
radiative efficiency rj = 0.2 were assumed in both cases. Taking 
fpeak = 300 s or 500 s has only minor effect on Tq in our case 
because the uncertainties in Ei so dominate. Rykoff et al. (2009) 
found Fo between ~100 and ~1000 in their ensemble of GRBs. 
The value of r in GRB 090726 thus falls into the lower half of 
this range and, given the uncertainty in its Ei so , it may even lie 
on the lower end. 

Ondfejov data show that the optical light curve achieves a 
power-law decay with a^o = -1.02 + 0.12 from t - To - 1730 s 
till the end of the series, which is a typical value for OAs. 
The X-ray observations of the afterglow of GRB 090726 started 
when the OA emission was already declining (see Fig.^d). 
The decline of the X-ray emission was power-law with a^x = 
-1.27 ± 0.09 before a break at t - Tq = 51+} 8 ks (Page et al. 
2009). Both ffdo and ardx can be regarded as identical at ~ 2cr 
level. The ratio of the optical and X-ray flux can be considered 
constant (Fig. [3]). All this suggests that both the optical and X- 
ray emission decayed simultaneously and that the spectral profile 
from X-ray to the optical band did not vary. This is true for both 
the time period before and after the break in the decaying light 
curve; no significant flares were observed in the data set of the 
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0.3-10 keV band. The SED in Fig. [4] suggests a break between 
the optical and X-ray band, with the X-ray part having a steeper 
slope than its extrapolation to the optical band, in the same sense 
as observed for OA of GRB 060418 by Molinari et al. (2007). 
Correction of the R band point by a very high Ar « 3.7 mag 
due to the extinction in the host galaxy would make it lie on the 
extrapolation of the X-ray spectrum in Fig.|4] However, the rela- 
tion between the optical extinction and X-ray absorption in the 
host is unknown. 

The break in the X-ray light curve occurred in the gap of the 
optical data, but the ratio of the optical and X-ray flux before 
and after the break shows that this event is consistent with an 
achromatic break. Following the approach of Ghirlanda et al. 
(2004), the Lorentz factor at the time of this break has already 
decreased to F^eak ~ 20. 

Although the intrinsic absorbing column of the X-ray after- 
glow of GRB 090726 is larger than the Galactic one (Page 2009), 
no changes of its H were observed during its power-law decline 
covering the part before the break in its light curve (Fig. [3]). This 
suggests a constant absorption which does not participate in the 
evolution of this afterglow. This can be explained if this absorp- 
tion comes from the regions more distant from this GRB, e.g. if 
this event is behind a molecular cloud. Any possible variations 
of the absorption caused by emission of this GRB or its after- 
glow must have occurred before the start of XRT observing, that 
is at t - To < 3000 s. The extinction of the optical emission 
depends on the dust content while the absorption of X-rays de- 
pends on the gaseous component of the absorbing medium. If 
only the dust is being destroyed during the early phase of the af- 
terglow or condensed during its later phase, then RX would vary, 
which is not observed in our case (Fig. [3]!. The constancy of both 
RX and H (within the observing scatter) suggests that neither the 
dust nor the gaseous component underwent any evolution during 
this OA event. 

We find that the power-law decay rate of OA of GRB 090726 
is comparable to that of the ensemble of OAs with homogeneous 
optical spectra at z = 0.36 - 3.5 (Simon et al. 2001, 2004). The 
^-corrected absolute magnitude (Simon et al. 2001) Mr m -22 
of OA of GRB 090726 at (t - r ) rest ~ 0.2 d places it below 
or at the lower end of the distribution of this ensemble of OAs. 
Mr of the ensemble span about 4 mag at this (f - 7o) re st- This is 
consistent with the synchrotron emission of the GRB afterglow 
fireball emission model by Sari et al. (1998) with very similar 
spectra but very different luminosity at a given (f- 7o) re st- In this 
phase, OA of GRB 090726 thus belongs to the least luminous 
ones and is similar to OA of X-ray flash XRF 030723 (Fynbo et 
al. 2004; Simon et al. 2006). 
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